8 Poincare sections 



The dynamical systems we study are of the form 




Systems of such equations describe a flow in phase space. 

The solution is often studied by considering the trajectories of such flows. 

But the phase trajectory is itself often difficult to determine, if for no other 
reason than that the dimensionality of the phase space is too large. 

Thus we seek a geometric depiction of the trajectories in a lower-dimensional 
space — in essence, a view of phase space without all the detail. 

8.1 Construction of Poincare sections 

Suppose we have a 3-D flow T. Instead of directly studying the flow in 3-D, 
consider, e.g., its intersection with a plane (x% = h): 




x 



1 



• Points of intersection correspond (in this case) to ±3 < on T. 



• Height h of plane S is chosen so that T continually crosses S. 



• The points Pq, Pi, P2 form the 2-D Poincare section. 



68 



The Poincare section is a continuous mapping T of the plane S into itself: 

Pk+i = T(P k ) = T [T(i\_i)] = T 2 (P fc _!) 
Since the flow is deterministic, Po determines Pi, Pi determines P2, etc. 

The Poincare section reduces a continuous flow to a discrete-time map- 
ping. However the time interval from point to point is not necessarily con- 
stant. 

We expect some geometric properties of the flow and the Poincare section to 
be the same: 

• Dissipation =^> areas in the Poincare section should contract. 

• If the flow has an attractor, we should see it in the Poincare section. 

Essentially the Poincare section provides a means to visualize an otherwise 
messy, possibly aperiodic, attractor. 

8.2 Types of Poincare sections 

As we did with power spectra, we classify three types of flows: periodic, 
quasiperiodic, and aperiodic. 

8.2.1 Periodic 

The flow is a closed orbit (e.g., a limit cycle): 

Po 





69 



Po is a fixed point of the Poincare map: 

p = r(p ) = r 2 (p ) = .... 

We proceed to analyze the stability of the fixed point. 

To first order, a Poincare map T can be described by a matrix M defined in 
the neighborhood of Pq: 



OXn 



J 3 

In this context, M is called a Floquet matrix. It describes how a point Po + 5 
moves after one intersection of the Poincare map. 



A Taylor expansion about the fixed point yields: 

dP 



dP 

P(P + 5) ~ T,(P ) + 



Since T(P ) = P , 



ox 2 



-S 2 , i = l,2 

Po 



T(P + 6) ~ P + M5 

Therefore 

t(t(Po + ^ - T(P + M5) 
~ T{P ) + M 2 5 
~ P + M 2 £ 

After m interations of the map, 

T m (P + 5)-P ~ M m 5. 
Stability therefore depends on the properties of M. 

Assume that 5 is an eigenvector of M. (There will always be a projection onto an 
eigenvector. ) Then 

M m 5 = \ m 5, 
where A is the corresponding eigenvalue. 

70 



Therefore 



|A| < 1 =4> linearly stable 
|A| > 1 =4> linearly unstable 

Conclusion: a periodic map is unstable if one of the eigenvalues of the 
Floquet matrix crosses the unit circle in the complex plane. 

8.2.2 Quasiperiodic flows 

Consider a 3-D flow with two fundamental frequencies f\ and f<z- The flow is 
a torus T 2 : 




The points of intersection of the flow with the plane S appear on a closed 
curve C. 

As with power spectra, the form of the resulting Poincare section depends on 
the ratio /1//2: 

• Irrational fi/f2- The frequencies are called incommensurate. The 
closed curve C appears continuous, e.g. 
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Xj 

— The trajectory on the torus T 2 never repeats itself exactly. 

— The curve is not traversed continuously, but rather 

T(C) = finite shift along C. 

• Rational /1//2. 

— fi and /2 are frequency locked. 

— There are finite number of intersections (points) along the curve 

— Trajectory repeats itself after n\ revolutions and rotations. 

— The Poincare section is periodic with 

period = m//i = n 2 // 2 

— The Poincare section contains just n\ points. Thus 

P l = T n \P l ) 

— Example, n\ = 5: 





8.2.3 Aperiodic flows 

Aperiodic flows may no longer lie on some reasonably simple curve. 
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In an extreme case, one has just a point cloud: 




This would be expected for statistical white noise. 

Deterministic aperiodic systems often display more order, however. In some 
cases they create mild departures from a simple curve, e.g. 




Such cases arise from strong dissipation (and the resulting contraction of 
areas in phase space). 

It then becomes useful to define a coordinate x that falls roughly along this 
curve, and to study the iterates of x. This is called a first return map. 

8.3 First-return maps 

First return maps are 1-D reductions of the kind of 2-D Poincare maps that 
we have been considering. 

Such maps are of the form 

Xk+i = f(x k ). 

We will study these extensively at the end of the course. 

We shall give particular attention to the following quadratic mapping of the 
unit interval onto itself: 

x k +i = 4:/ix k (l - Xk), </i<1. 
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The mapping is easily described graphically. The quadratic rises from x = 0, 
falls to x = 1, and has its maximum at x = 1/2, where it rises to height fi. 



Consider, for example, the case fi = 0.7: 



fix) 



1 




f(x Q ) =Xj 



x, 
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X 



1 



X 



Eventually the interations converge to x = x, which is where the diagonal 
(the identity map Xk+i = Xk) intersects f(x). 

Thus x is a fixed point of /, i.e., 



Another fixed point is x = 0, since f(0) = 0. 

However we can see graphically that x = is unstable; iterates initiated at 
xq = e still converges to x. 

Thus x = is an unstable fixed point, while x = x is stable. 
What determines stability? Consider graphically the case fi = 0.9: 



x = f(x) 



74 



fix) k 
A 

1 ~ 




1 x 



We infer that the slope f'(x) determines whether x is stable. We proceed to 
show this formally. 

Suppose x* is any fixed point such that 

x* = /(**). 

Define 

Xk = x* + Ek, £k small. 
In general, our mappings are described by 

Xk+i = f(x k ). 

Then 

X*+£ k+1 = f{x* + £ k ) 

= f(x*) + f'(x*) £k + 0(e 2 k ) 

Therefore 

£k+i ^ f'{x*)e k . 

Thus 

\f'(x*)\ < 1 stability. 
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8.4 1-D flows 



It is instructive to compare the stability of 1-D maps to the stability of the 
continuous 1-D flow 

x = f(x). 

Consider, for example, a particular f(x) that looks like 

x 





yf(x) 


V ^ 


- x 


^^^^ 


^, — x 2 







Clearly x\ and X2 are fixed points. Which are stable? 

The arrows show the direction of flow on the x-axis, i.e., the sign of x = f(x). 
Thus X\ is stable while x^ is not. 

Stability of a fixed point x* is therefore determined as follows: 

f\x*) < => stable 

f'(x*) > => unstable 

Whereas the stability of a continuous 1-D flow / depends on the sign of /', 
the stability of a 1-D map depends on the magnitude \f'\. 

In higher dimensions this same distinction holds for the eigenvalues A of 
the Jacobian (which, in the case of mappings, we have called the Floquet 
matrix). That is, the sign of Re (A) determines the stability of flows, whereas 
the magnitude |A| is the relevant quantity for maps. 
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8.5 Relation of flows to maps 



We now consider explicitly how flows may be related to maps. 



8.5.1 Example 1: the van der Pol equation 



Consider again the van der Pol equation 

d 2 9 _ .,d0 



dt2+£{e >-i )ft+ e = 

Recall that for e > the rest position is unstable and that the system has a 
limit cycle. 

We draw a ray emanating from the origin, and consider two representative 
trajectories initiating and terminating on it: 



1 
























Let Xk be the position of the kih intersection of the trajectory with the ray. 
There is then some mapping / such that 

Xk+i = f(x k ). 
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The precise form of f(x) is unknown, but physical and mathematical reason- 
ing allows us to state some of its properties: 



• / maps Xk to a unique Xk+i- 

• / is continuous. 

• f'(x) > 1 near the origin (divergent spirals). 

• f'{x) < 1 far from the origin (convergent spirals). 

• f'(x) > for all x > (since f(x + 5) > f(x)). 



The simplest form of / is therefore a curve rising steeply from the origin, 
followed by a gentle upward slope: 




X k+1 ~ X k 



X 



By continuity, there must be a stable fixed point x* characterized by 

x* = f( x *) and f'(x*) < 1. 
Thus x* gives the effective radius of the stable limit cycle. 



8.5.2 Example 2: Rossler attractor 



Consider the following 3-D flow (the Rossler attractor): 

x = —y — z 
y = x + ay 
i = b J rz(x — c) 
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a, 6, and c are fixed parameters. 
Numerical solutions yield the time series x(t): 

Figure IV. 11, BPV 

The time series z(t): 

Figure IV. 11, BPV 
A Poincare section in the x-y plane: 

Figure IV. 12a, BPV 

And a 3-D perspective of the flow: 

Figure 12.1, unknown source 

The time series display great irregularity, but the Poincare section and the 
full flow display some order. 

Consider now another Poincare section, in the plane 

y + z = 0. 

From the Rossler equations, we identify this plane with extrema in the time 
series x(t), i.e., each intersection of the plane corresponds to 

x = 0. 

Consider a sequence Xk of such extrema, but only when the extremum is a 
maximum (local peak) of x(t). 

Then plot Xk+i vs. x^. 

Figure IV. 10, BPV 

Conclusions: 

• The 1-D map, even more so than the Poincare section, reveals that the 
flow contains much order. 

• The time series, however, displays no apparent regularity. 
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This is the essence of deterministic chaos. 

We proceed to show how such Poincare sections and 1-D maps can be con- 
structed from experimental data. 

8.5.3 Example 3: Reconstruction of phase space from experimental data 

Suppose we measure some signal x(t) (e.g., the weather, the stock market, 
etc.) 

In most cases it is unlikely that we can specify the equations of motion of the 
dynamical system that is generating x(t). 

How, then, may we visualize the system's phase space and its attractor? 

The (heuristic but highly successful) idea, originally due to Santa Cruz un- 
dergraduates Packer, Crutchfield, Farmer, and Shaw (1980), is to measure 
any 3 independent quantities from x(t). 

For example: 

• x(t), x(t + r), x(£ + 2r); t large enough for "independence," i.e., beyond 
an autocorrelation time. This is the most popular; it is known as the 
method of delays. 

• x(t), x(t), x(t) (where the derivatives are finite differences Xk — Xk-i, 
etc.). 

Such a representation of the attractor is not identical to the "real" phase 
space, but it should retain similar geometric properties. 

We use the Rossler attractor as an example. 

• Comparison of the projection of the trajectories in the x-y plane to the 
projection of the trajectories in the x-x plane: 

Figure IV. 12, BPV 



80 



• Example of different values for r using the method of delays: 



Figure 78, Schuster 



With the method of delays, r is typically the period of the forcing, or the 
period of a characteristic limit cycle. 

Here we have discussed only qualitative, geometric properties. We shall see 
that the various representations also yield similar quantitative properties 
(e.g., measures of Lyaponov exponents). 
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